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P_{ , Abstract 



We study the existence of a unique stationary distribution and ergodicity for a 2-dimensional 
affine process. The first coordinate is supposed to be a so-called a-root process with a € (1,2]. 
The existence of a unique stationary distribution for the affine process is proved in case of a G (1,2]; 
further, in case of a = 2, the ergodicity is also shown. 



rn ' 1 Introduction 
in 
(N 

' We consider the following 2-dimensional affine process (affine two factor model) 

o ■ 

cn; {dYt = {a-bYt)dt+ ^dLt, t^O, 



(1.1) 



dXt = {m-eXt)dt + ^/YtdBt, t ^ 0, 



^ , where a > 0, 6, 6',?7i G M, a G (l5 2], {Lt)t^o is a spectrally positive a-stable Levy process with 
I Levy measure Ca2~^~"l|^>o} with Ca '■= {aT{—a))~^ (where F denotes the Gamma function) in 
case aE(l,2), a standard Wiener process in case a = 2, and {Bt)t^o is an independent standard 
Wiener process. Note that in case of a = 2, due to the almost sure continuity of the sample paths 
of a standard Wiener process, instead of y/Yi~ one can write \/Yt in the first SDE of (|l.ip . and 
Y is the so-called Cox-Ingersol-Ross (CIR) process; while in case of a S (1,2), Y is called the 
a-root process. Note also that the process (It)f^o given by the first SDE of (jl.ip is a continuous 
state branching process with immigration with branching mechanism bz + ^z"^, z ^ 0, and with 
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immigration mechanism az, z ^ (for more details, see the proof of Theorem 13.11 part (i)). Chen 
and JosHn [7] have found several applications of the model (jl.ip with a = 2 in financial mathematics, 
see their equations (25) and (26). 

The process (Y, X) given by (jl.ip is a special affine process. The set of afRne processes contains a 
large class of important Markov processes such as continuous state branching processes and Orstein- 
Uhlenbeck processes. Further, a lot of models in financial mathematics are affine such as the Heston 
model |16] . the model of Barndorff-Nielsen and Shephard [3] or the model due to Carr and Wu [6]. A 
precise mathematical formulation and a complete characterization of regular affine processes are due 
to Duffie et al. [10]. Later several authors have contributed to the theory of general affine processes: to 
name a few, Andersen and Piterbarg [1], Dawson and Li [9], Filipovic and Mayerhofer |11] . Glasserman 
and Kim [H], Jena et al. [E] and Keller- Ressel et al. [23] . 

This article is devoted to study the existence of a unique stationary distribution and ergodicity of 
the affine process given by the SDE (jl.ip . These kinds of results are important on their own rights, 
further they can be used for studying parameter estimation for the given model. For the existing 
results on ergodicity of affine processes, see the beginning of Section [3l 

Next we give a brief overview of the structure of the paper. Section [2] is devoted to a preliminary 
discussion of the existence and uniqueness of a strong solution of the SDE (jl.ip by proving also that 
this solution is indeed an affine process, see. Theorem 12. 2[ In Section [3] we prove the existence of a 
unique stationary distribution for the affine process given by (jl.ip in both cases a G (1, 2) and a = 2, 
provided that a > 0, b > and ^ > 0, see, Theorem 13. 1[ In Section U in case of a = 2, we prove 
ergodicity of the process in question provided that a > 0, 6 > and ^ > 0, and we also show that 
the unique stationary distribution of the process is absolutely continuous, has finite (mixed) moments 
of any order by calculating some moments explicitly, too, see Theorems 14.11 and 14. 2j respectively. 

In a forthcoming paper we will use our results for studying parameter estimation for the process 
given by the SDE (fTT]) . 



2 The affine two factor model 

Let N, M and M-|_ denote the sets of positive integers, non-negative integers, real numbers and 
non- negative real numbers, respectively. By and ||^|| we denote the Euclidean norm of a vector 
X G M™' and the induced matrix norm ||^|| = sup{||^a;|| : x € M™, ||x|| = 1} of a matrix A G M"^™, 
respectively. By C'^{R+ x M,M), C^{R+ x M,M) and C^{R+ x M,M), we denote the set of twice 
continuously differentiable real- valued functions on xM, the set of twice continuously differentiable 
real-valued functions on x M with compact support and the set of infinitely differentiable real- 
valued functions on IR+ x R with compact support, respectively. Convergence in distribution will 
denoted by — >. 

Let {^,T,{J^t)t^o,^) be a filtered probability space satisfying the usual conditions, i.e., {^,T,F) 
is complete, the filtration {J-t)t^o is right-continuous and J-q contains all the P-null sets in J^. Let 
{Btjt^o be a standard (-Fi)i^o-Wiener process and {Lt)t^o be a spectrally positive (J"4)fj.o-stable 
process with index as(l,2]. We assume that B and L are independent. If a = 2, we understand 
that L is a standard (/"t)(j>o-Wiener process. If aE(l,2), we understand that L is a (J-()t^o-L6vy 
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process with Levy-Khintchine formula 



E(e^"^i) = exp 1^ (e*"^ - 1 - iuz)C a^z'^''^ dzj , u& M, 

where Cq, = {aT{—a))~^ . Recah that in case of a G (1,2) the Levy-Ito representation of L takes 
the form 

Lt= I I ziV(ds,dz), t ^ 0, 

J{Q±\ J(0,oo) 



{0,t] J(0,oo) 

where A^(d,s, dz) is a compensated Poisson random measure on (0, oo)^ with intensity measure 
CcZ-^-"l|2>0} dsdz 



2.1 Remark. We shed some Hght on the definition of the stochastic integral with respect to the 
spectrally positive a-stable process L in the first SDE of (|1.1|) in case of a € (1,2). By Jacod 
and Shiryaev |18|. Corollary II. 4. 19], L is a semimartingale so that Theorems 1.4.31 and 1.4.40 in 
Jacod and Shiryaev |18] describe the classes of processes which are integrable with respect to L. A 
more accessible integrability criteria is due to Kallenberg \20\ Theorem 3.1]. Roughly speaking, a 
predictable process V is locally integrable with respect to L (i.e., the stochastic integral J^^VgdLs 
exists for all t ^ 0) if and only if Jq \ Vs\°' ds < oo almost surely for all t ^ 0. For the construction 
of stochastic integrals with respect to symmetric a-stable processes, see also Rosinski and Woyczynski 
[3H Theorem 2.1]. Another possible way is to consider the stochastic integral with respect to L 
as a stochastic integral with respect to a certain compensated Poisson random measure, see the last 
equality on page 230 in Li [24]. 



The next proposition is about the existence and uniqueness of a strong solution of the SDE (jl.ip . 

2.2 Theorem. Let (lo,^o) be any J^Q-measurable random vector such that P(lo ^ 0) = 1. Then 
for all a G ffi+, b,m,9 £ M, a £ (1,2], there is a (pathwise) unique strong solution (Yt, Xt)t^o of 
the SDE ([□]) such that F{Yt ^ 0, V t ^ 0) = 1. Further, we have 

(2.1) = e-^(*-^) (Y, + aj\-'^'-^^du + j\-''^'-''^0^dLu^ 
for ^ s ^ t, and 

(2.2) A:t = e-^(*-") (^X, + m^*e-^("-"M'u + ^*e-^(^-")v^dB„^ 

for ^ s ^ t. Moreover, (Yt, Xt)t^o is an affine process with infinitesimal generator 
i-^f)iy,x) = (a - by)f[{y,x) + {m - ex)f2{y,x) + ^yf2,2iy,x) 

(2.3) 



in case of a £ (1,2), and 



(2.4) {Af){y,x) = {a- by)f[{y,x) + (m - 9x)f^{y,x) + \y{fUy,x) + fl2{y,x)) 

in case of a = 2, where {y,x) G R+ x M, / E Cc(M+ x R,R), and f-, i = 1,2, and fl'j, 
i,j G {1,2}, denote the first and second order partial derivatives of f with respect to its i-th and 
i-th and j-th variables. 
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Proof. By Theorem 6.2 or Corollary 6.3 in Fu and Li ^13j (for the case a E (1, 2)) and by Yamada and 
Watanabe theorem (see, e.g., Karatzas and Shreve |21l Proposition 5.2.13]) and Ikeda and Watanabe 
[171 Example 8.2, page 221] (for the case a = 2), there is a (pathwise) unique non-negative strong 
solution (Yt)t^o of the first equation in (11. ip with the given initial value Yq. An application of the 
Ito's formula to the process (e**Yf)f^o shows that 

d(e**yj) = be''^Yt dt + e''^dYt = be^% dt + e** ((a - bYt) dt + ^/Y^dLt) 

= ae^^ dt + e'"' ^/Ytl dLt, t ^ 0. 



Hence, if ^ s ^ then 



% = e'^Ys + a f e^" dn + /* e^" 

J s J s 



and (I2JD follows. 

Let {Xt)t>o be defined by ()2.2p . By Ito's formula, we have 

dXt = -ee-^*(^Xo + m e^'ds + e^* y^dS,) dt + e-^\me^^ dt + e^*v^d5t) 

= -ee-^*(^Xo + m e^'ds + e^'y^sdBs^dt + mdt + ^/YtdBt, t ^ 0. 

Then the process {Xt)t>o is indeed a strong solution of the second equation in (jl.ip . giving the 
existence of the strong solution of (jl.ip . Conversely, if {Xt)t^o is a strong solution of the second 
equation in (jl.ip . by an application of Ito's formula to the process {e^^Xt)t^o, we have 

d{e^^Xt) = Oe^^Xt dt + e^'MXj = Oe^^Xt dt + e^*((m - OXt) dt + ^/YtdBt) 

= me^*dt + e^*V^dSt, t ^ 0. 

If ^ s ^ t, then 

e^*Xf = e^'X, + m [ e^" dn + / e^^^l^d^^. 



and hence {Xt)t>o must be given by (j2.2p . This proves the pathwise uniqueness for the second 
equation in (jl.ip . and hence the pathwise uniqueness for (jl.ip . 

Now we turn to check that (It, Xt)t^o is an affine process with the given infinitesimal generator. 
In case of a = 2, by Ito's formula, for all / G (M+ x ]R,M) we have 

fiYt,Xt) = f{Yo,Xo)+ f f[(Ys,Xs)VYsdLs+ T /2(^., V^di?, 

Jo Jo 

+ f f[(Ys,Xs){a-bY,)ds+ f f^{Ys,Xs){m-eXs)ds 
Jo Jo 

+ \ [J^ fi,i{Ys, Xs)Ys ds + 1^ fl^{Ys,Xs)Ys d^ 

= f{Yo,Xo)+ [\Af){Ys,Xs)ds + Mtif), t^O, 
Jo 
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where 



Mt{f):= f f[{Ys,Xs)^/YsdLs+ T /^(n, X,)^/^ dS„ t ^ 0, 
Jo Jo 



and Af is given by (j2.4p . If E(yo) < oo, then {Mt{f))t^o is a square integrable martingale with 
respect to the filtration {Tt)t^o, since 

/ E{{f[{Ys,Xs))^Ys)ds!^Ci [ E(n)ds<oo, t ^ 0, 
Jo Jo 

[ E{{f^{Ys,X,))^Ys)ds^C2 [ E(n)ds<oo, t ^ 0, 
Jo Jo 

with some constants Ci > and C2 > 0, where the finiteness of the integrals follow by Proposition 
3.2 in Barczy et al. [2]. 

In case of a G (1, 2), by the Levy-Ito representation of L, we can rewrite the SDE (|l.ip into the 
integral form 

(Yt = Yo + f^{a - bY,) ds + z ^/YZN(ds, dz), t ^ 0, 

\Xt = Xo + /d(m - 9Xs) ds + j/^VZdBs, t ^ 0. 

By Ito's formula, for ah / G C^(M+ x R,R) we have 

f{Yt,Xt) = f{Yo,Xo)+ f f[{Ys,Xs){a-hY,)ds+ f /^{Y^, Xs){m - OX^) ds 

Jo Jo 

+ 1 f f^2iys,Xs)Ysds+ f f^{Ys,Xs)VYsdBs 

^ Jo Jo 
+ r r (fi^s- + z VlT, X,J) - f(Ys^, Xs^)) N{ds, dz) 



10 JO 

ft roo 



^0 



+ / / (f{Y, + z^Ys,Xs)-fiYs,X,)-z^Ysf[iYs,Xs))Co,z~'~''dsdz 



where 



= /(yo,^o)+ [\Af){Ys,X,)ds + Mt{f), t^O, 
Jo 

Mtif) := f f^{Ys,Xs)VYsdBs 
Jo 



+ [f{Ys- + z ^1^, Xs-) - f{Ys-,Xs-)) N{ds, dz), t ^ 0, 

and, by the change of variable z ^ := z, 

{Af){y,x) := {a - by)f[{y,x) + (m - ex)f^{y,x) + ^yfi^iy,^) 



+ y [ I" {f{y + z,x)-f{y,x)-zf[{y,x))CaZ-'-''dsdJ 
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for {y,x) eR+xR and f e C^{R+ xR,R). Next we check that if E(yo) < oo, then {Mt{f))t^o is 
a martingale with respect to the filtration {Tt)t^o- The martingale propery of Jg f2{Ys,Xs)\/YsdBs, 
t ^ 0, follows in the very same way as in the case of a = 2 using that there exists some constant 
C3 > such that E{Yt) ^ C3(l + E(yo)e"''*/°) for ah t ^ 0, see, e.g., formula (2.7) or Proposition 
2.10 in Li and Ma [25]. Now we turn to check that 

M^if) := 1^ (^f{Ys^ + z 0^,Xs^) - f{Ys^,Xs^)) N{ds,dz), t ^ 0, 

is a martingale. For all n E N, let 

'"(/) := (/((n- A n) + z A n, - /(y,_ A n, iV(ds, dz), 

Mj''"(/) (/((F.- A n) + 2 ^y,_ A n, - /(y,_ A n, iV(ds,dz). 

By Taylor's theorem, there exists some random variable Cn (depending also on s and z, but not 
denoting this dependence) such that 

/((n_ A n) + z V>^r^, X,_) - /(y,„ A n, X,_) 

= /( ((y,_ A n) + C„z \/ys- A n, z v^y,_ A n. 

Since 

(/((y,_ An) + z^/l^r^,X,_) -/(y,_ An,X,_))c7„z-i-"dsdz 
= yi y E ((y, A n) + Cn^ \/y A n, X,) ^y, A n) C^z'^ds dz 

/■oo J. 

^ dn^l'^t / dz = dn^/'^ < 00, M+, 

Jl Oi-l 

with some constant C4 > 0, by Lemma 3.1 in Chapter I and page 62 in Ikeda and Watanabe [17] . 
we get (M^^'"(/))t^o is a martingale with respect to the filtration (-Fj)t^o- Further, since 

^(i) i) (^((^-- ^ ^) + 2 V^'.- A n, X,_) - /(y,_ A n, X,„)) ' C^z-^^^ds dz 

^ C5 r E((y, A n)2/°) ds /"^ zi-° dz < C^n^/'^-L_ < 00, i e M+, 
Jo Jo 2 - a 

with some constant C5 > 0, by pages 62 and 63 in Ikeda and Watanabe [T7]. we get {Mf'^{f))t^Q is a 
square integrable martingale with respect to the filtration {J-'t)t^o- This yields the martingale property 
of {Ml{f))t^o with respect to the filtration {Ft)t^Q. Indeed, for all nGN and / G C^(R+ x M, R), 
let 

%"(/) := MlU) - Mf'"(/) - Mf'"(/), t G R+. 
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Then, by Taylor's theorem, 

ct roc 



ritU) = 11" l{y,_>„}(/(ys- + z v^,X,_) - f{Ys-,Xs-)) N{ds,dz) 

f't f'OO ^ 

Jo Jo 

rt roo 

= / / l{Y^_>n}f[{Ys- + Cz^/Y:I,Xs-)^/Y;ZdL,, tGR+ 
Jo Jo 

with some random variable Using Lemma 2.8 in Li and Ma [25], there exist some constants > 
and C7 > such that for all t G ]R_|_, 



E ( sup CeE n^*l|y^>„|(/{(n + C2V^,X,))"y,ds) 



=^C7^*E(i{y^>„|y,)ds) ^ , 

where the last inequality follows by Jensen inequality using also that /' is bounded. If E(yo) < 00, 
then there exists some constant Cg > such that E(lt) < 6*8(1 + E(yo)e~^*/"), t G M+, see, e.g., 
formula (2.7) or Proposition 2.10 in Li and Ma |25j . Hence 

Te (l|y^>„}y,) ds ^ Cs Al + E(>o)e-'*/")ds < 00, n e N, 
Jo Jo 

and then, by dominated convergence theorem, 



(2.5) lim E ( sup |?7^(/)| =0, i G 

.sG[0,t] I 



n— >oo 



This yields that (Af/(/))t^o is a martingale with respect to the filtration {Ttjt^o- Indeed, for all 
^ s s$ t and A G 7"^, by ([23]) . we have 



E(|r?r(/)lA|)^E(|r?r(/)|)^0 as n ^ 00, 
E(|7?«(/)1^|)^E(|<(/)|)^0 as n^oo. 

Hence 

lim E((m2'"(/) + M,3'"(/))1^) =E(Mi(/)l^), 

n— >oo 

(2.6) 

lim E((m2'"(/) + m3'"(/))1^) = E(Ml(/)l^). 

n— >oo 

Since (M^^'"(/) + M^'^ {f))t^Q is a martingale with respect to the filtration (-Ft)j5.o for all n G N, 
we have 

E((M,2'"(/) + Mi^'"(/))l^) =E((m2'"(/) + m3'"(/))1^), nGN, 
which, by ()2.6p . yields that 

E(m1(/)1^) =E(A./i(/)lA) 
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for all O^s^t and A Tg, i.e., (M^^(/))(j>o is a martingale with respect to the filtration 



Finally, we check that the transition semigroup {Pt)t^o with state space M+ x M correspond- 
ing to {Yi;,Xt)t^o is a regular affine semigroup having infinitesimal generator given by (j2.3p and 
(|2.4p according to the cases a G (1)2) and a = 2. With the notations of Dawson and Li [9], 
(0, {aij)lj^^, {hifi^i, {Pi,j)lj=i, 0, /u) and (0, {aij)lj^;^, (6i)i=i' (Aj)?j=i> 0. 0) are sets of admissible 
parameters according to the cases a G (1,2) and a = 2, where 



,"2,1 "2,2, 




[a,m) 



for a G (1,2], 



/3i,i /3i,2 




for a G (1,2], 



n{Ay,dx) := CaV ^ "dy5o(dx), 

where 5q denotes the Dirac measure concentrated on G 
6.1 in Dawson and Li [9] holds, since 



Indeed, condition (vi) of Definition 



(|y| Ay2)^(dy,dx) + 



{\x\ A x'^)fi{dy, dx) 



{\y\/\y')y-'-''dy6o{dx) = a 



POO 

Jo 



^y'^)y 



-l-a 



dy 



y 



' dy < oo. 

/o Ji 

Hence Theorem 2.7 in Duffie et al. |10j (see also Theorem 6.1 in Dawson and Li [9|) yields that for 
these sets of admissible parameters, there exists a regular afhne semigroup {Qtjt^o with infinitesimal 
generator given by (12. 3p and (12. 4p according to the cases a G (1,2) and a = 2. By Theorem 2.7 
in Dufhe et al. [1^, C^(M+ x R, M) is a core of the infinitesimal generator corresponding to the 
affine semigroup {Qt)t^o- Since we have checked that the infinitesimal generators corresponding to 
the transition semigroups {Pt)t^o and {Qt)t^o (defined on the Banach space of bounded functions 



on 



coincide on 



+ X M, M), by the definition of a core, we get they coincide on the 
Banach space of bounded functions on IR+ x R. This yields that (It, Xt)t^o is an affine process with 
infinitesimal generator ()2.3p and ()2.4p according to the cases a G (1, 2) and a = 2. We also note 
that we could have used Lemma 10.2 in Duffie et al. jlO] for concluding that {Yt,Xt)t^o is an affine 
process with infinitesimal generator (j2.3p and ()2.4p according to the cases a G (1,2) and q = 2, 
since we have checked that if E(yo) < oo, then {Mt{f))t^o is a martingale with respect to the 



filtration (J^t)t>o for any f £ 



in both cases a G (1,2) and a = 2. 



□ 



2.3 Remark, (i) Note that in Theorem 12.21 the assumption that {Yq,Xo) is J^o-™6asurable yields 
that {Yo,Xq) is independent (of the increments) of {Lt,Bt)t^o- 

(ii) Note that in Theorem 12.21 it is the assumption a G M+ which ensures P(lt ^0, V t ^ 0) = 1. 
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3 Stationarity 



The study of existence of stationary distributions for affine processes in general is currently under 
active research. 

In the special case of continuous state branching processes with immigration the question of exis- 
tence of a unique stationary distribution has been well-studied, see Li [241 Theorem 3.20 and Corollary 
3.21] or Keller-Ressel and Mijatovic [22\ Theorem 2.6]. 

Glasserman and Kim |14t Theorem 2.4] proved existence of a unique stationary distribution for 
the process 



(3.1) 



dYt = ia-bYt)dt + VYtdLt, t^O, 



[dXt = -eXt dt + y/l + cjYtdBt, t ^ 0, 

where a>0, & >0, 0>O, u^O and L and B are independent standard Wiener processes. 

The following result states the existence of a unique stationary distribution of the affine process 
given by the SDE (jl.ip for both cases a S (1,2) and a = 2. 

3.1 Theorem. Let us consider the 2-dimensional affine model (11. ip with a > 0, 6 > 0, m G M, 
6 > 0, and with a (random) T^-measurahle initial value {Yq,Xq) such that P(lo ^ 0) = 1. Then 

(i) {Yf,Xt) — > {Y^,X^) as t — t- oo, and the distribution of (^005^00) is given by 
(3.2) E (e-AiFoo+iA^x^^ =exp|-a^"^;,(Ai,A2)ds + «^A2| 

for (Ai,A2) G K+ X R, where Vt[\i,\2), t ^ 0, is the unique non-negative solution of the 
(deterministic) differential equation 



(3.3) 



'|^(Ai, A2) = -bvtiXi,X2) - ^ivt{Xi,X2)r + ke~^''Xl t ^ 0, 
vo{Xi, A2) = Ai. 



(ii) supposing that the random initial value {Yq,Xq) has the same distribution as (yoo-'^oo) given 
in part (i), we have {Yf, Xt)t^o is strictly stationary. 

Proof, (i): Step 1. In this step we give some representations of the affine transition semigroup 
(-Pt)i>o with state space M+ x M corresponding to the process given by the SDE (jl.ip . By Theorem 
6.1 in Dawson and Li [S] and Theorem 12.21 the transition semigroup {Pt)t^o is given by 

roo roo 

(3.4) / / e<"'«>Pj((yo,xo),dO =e«^0'^o)'^*(")>+'^*(") 

Jo J -00 

for ueU, {yo,xo) G E+ x M, t^O, where [/ := C_ x (iR) with 

C_ := {zi + iz2 : zi E (-00, 0], Z2 £ M}, iR := {iz2 : Z2 £ R}, 
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and for all u = {ui,U2) G U, we have iptiu) = {i't^Huj^e ^*W2)> i ^ 0, where ^pl^\u), t ^ 0, is a 
solution of the generalized Riccati equation 



(3.5) 
and 



^^p[ 



i^oHu) = Ui, 



Jo 



'0 



Here for a G (1,2], the (complex valued) functions F and R are given by 

F{u) = aui + mu2, R{u) = —bui -\ 1 — -, u={ui,U2)€U. 

a 2 

Indeed, in case of a G (1,2), the formula for R{u), u £U, can be checked as follows. By Dawson 
and Li [9], in case of a G (1, 2), 



u 



2 /"oo /-oo 
2 



J~OD 



„,2 foo 

= -bui + ^ + c« / (e"^«i - 1 - ^zieoer'-" 

^ Jo 

(-til)" 

a 2 

where for the last equality we used that l/T{—a) = a{a — 1) /r(2 — a), the imaginary part of —iui^i 
is non-negative and that 



r(2 - a) 70 

for all complex numbers z with non-negative imaginary part, see, e.g., Zolotarev \37\ pages 67 and 
68]. 

Note also that for all u = {ui,U2) G U and t ^ 0, the real part of ipl^\u) is less than or equal 
to (compare also with Remark 2.2 in Duffie et al. |10j). and, in addition, if ui G M such that 
ui ^ 0, then Tpl^\u) G M with ^ 0. Moreover, for all t ^ 0, we have 

(ptiu) = / {aip^^\u) + me~^^U2)ds = a / ipl^^u) ds + mu2 ^ ■ 



In fact, one can give a simplified characterization of the transition semigroup {Pt)t^o by 

(3.6) e-^'^'+'^^^^Pt{{yo, xq), dO = exp {-yo^;t(Ai, A2) + ixoe'^^Aa + gt{Xi, X2)} 
for (Ai,A2), {yo,xo) G M+ x M, where 

(3.7) 5't(Ai,A2):= / {-aVs{Xi, X2) + irne~^^ X2) ds = -a Vs{Xi, X2) ds + imX2 



1 - e-^* 
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^(Ai, A2) = -bvtiXi, A2) - ^{vt{Xi,X2)r + ^e-2^*Ai, t ^ 0, 
■yo(Ai, A2) = Ai 



and ft(Ai, A2), t ^ 0, is the unique non-negative solution of the differential equation 
(3.8) 

in case a G (1,2]. Indeed, by (j3.4p with ui := — Ai and U2 := iA2, we have 

/"OO /"OO 

Jo J-oo 

= exp jyoV't (-Ai,iA2) + i2;oe~^*A2 + a y Vi^^ (-Ai, iA2) ds + zmA2 ^ 

where 

^(-Ai,a2) = -^V'f ^(-Ai,iA2) + -(-Vf ^(-Ai,a2))" + -e-2^*(iA2)^ t ^ 0, 

with V'o^'' (~Ai, iA2) = — Ai incase a G (1,2]. Introducing the notation ft(Ai,A2) := —ilj^^\—Xi,iX2), 
we have the differential equation (|3.8p for a G (1, 2]. Note also that 

vt{Xi, A2) ^ for all (Ai, A2) e M+ x M, 

since il)^\u) ^ for ui ^ 0. The proof of the uniqueness of the solutions of the differential equation 
(Cauchy problem) (|3.8|) is based on the global version of the Picard-Lindelof existence and uniqueness 
theorem: if / : D — )• R is a continuous function on a connected, open set L) C satisfying the 
global Lipschitz condition in its second variable 

\f{t, zi) - fit, Z2)\ < C\zi - Z2I, [t, zi), (t, Z2) G D, 

with some constant C > 0, then for all (to,^o) £ D, the Cauchy problem 

z'{t) = f{t,z{t)) with initial value z{tQ) = zq 

has a unique solution z defined on a maximal interval of the form (t_,t+), where t_ < 
t- G [—00,00), t+ G (—00,00], and (t,z{t)) leaves any compact subset of D as t 4, t_ and 
i t ^+) see, e.g., Hartman [151 Chapter II, Theorems 1.1 and 3.1] and Walter [361 Chapter 3]. Given 
(Ai, A2) G M+ X M, for all n ^ Ai, n G N, let Dn := (-l,oo) x {-n,n) and let /„ : -D„ M be 
defined by 

fn{t, z) := -bz - -z" + ^e-2^*A^, [t, z) G with a G (1, 2]. 

Further, let Iq := and zq := Ai. Then Dn is open and connected, /„ is continuous and 
satisfies the global Lipschitz condition in its second variable since, by the mean value theorem, for all 
(ii,^i), (^2,2:2) G -D we get 

|/„(tl,Zl) - fn{t2,Z2)\ ^ b\zi - Z2I + -\zf - ^2 I b\zi - Z2I + - Z2\. 

a 

By the global version of the Picard-Lindelof existence and uniqueness theorem, there is a unique 
solution Zn of the Cauchy problem 

^n(*) = Znit)) with initial value z„(to) = zq, 
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defined on a maximal interval of the form + where (in)- < (^n)+> (^n)- £ [—00,00) 

and (t„)+ E (— 00, +00]. Further, the solution Zn leaves every compact subset of D„ which implies 
that {tn)+ = 00 for all n G N. This yields the uniqueness of the solutions of the differential equation 
(Cauchy problem) (|3.8p . 

Step 2. We show that 
(3.9) t;t(Ai,A2) 5SM(Ai,A2)(l + i)max(e"^^*,e"^*), * ^ 0, (Ai, A2) G M+ x M, 

where 

A? 



M(Ai,A2) : = 



Ai + ^ ifM20. 
Ai + ^ if 6 = 26. 



2 

The proof is based on the following version of comparison theorem (see, e.g., Volkmann [35] or Lemma 
B.3. in Filipovic et al. |12j): if S : M-|_ x M — ^ M is a continuous function which is locally Lipschitz 
continuous in its second variable and p, q : M+ — )• M are differentiable functions satisfying 

p'it) ^ Sit,p{t)), t ^ 0, 
p{0) q{0), 

then p{t) ^ q(t) for all t ^ 0. Using this one can check that i'j(Ai,A2) ^ uj(Ai,A2) for all t ^ 
and (Al, A2) € x M, where for all (Ai, A2) G M+ x M, ut{\i, A2), t ^ 0, is the unique solution of 
the differential equation 



(3.10) 



||^(Ai,A2) = -fent(Ai,A2) + ie-2eiA2, t ^ 0, 



['Uo(Ai, A2) = Al. 
Further, one can verify that 



(3.11) ut{Xi,X2) 



\_, I _^2__] „~bt ^2 ^~20t :f i, no 

^ 2{-b+29) J 2{-b+29)^ LI ^ Zt7, 

Ai + ^t)e-^* if b = 26. 



Indeed, the general solution of the homogeneous differential equation ^(Ai,A2) = — 6nt(Ai,A2), 



dt 

-bt 



t ^ 0, takes the form nt(Ai,A2) = Ce , t ^ 0, where C S M, and it can be checked that a 
particular solution of the inhomogeneous differential equation (j3.10p (without the initial condition) is 



iif(Ai, A2 



2(-6+26»)^ n o ^ zc, 

^ie-2^* if b = 2e. 



Hence the general solution of the differential equation (|3.1Up (without the initial condition) takes the 
form 

jLme-^'* if b + 26, 



Then taking into account the initial condition tio(Ai, A2) = Ai, we have (j3.1ip . 



fte-2^* if b = 26. 
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Finally, by (f37TT]) . we readily have (f3^ in case of & = 20. If 6 > 26*, then 
and if <b <26, then 
as desired. 



2{-b + 2e) J \ 2{-b + 2e) 



Step 3. By the continuity theorem and p.6p . to prove (i), it is enough to check that for all 
(Ai,A2), iyo,xo) G M+ X M, 

lim {-yoVtiXi, A2) + ixoe~'^^X2 + gt{Xi, X2)) = -a Vs{Xi,X2) ds + i— A2 
(3.12) -^0 

='■ 9oo{Xi, A2), 

and that the function x R 9 (Ai,A2) 1— >• 5oo(Ai,A2) is continuous. Indeed, using (|3.6p and the 
independence of (Yoi-'^o) and {L,B), the law of total expectation yields that 

i/ — " 00 

expl-yoVt{Xi,X2) + ixoe~^^X2 +5t(Ai,A2)[ P(yo,Xo)(dyo, dxo) 

for all (Ai,A2) G IR+ x M, where P(yg Xo) denotes the distribution of {Yq,Xq) on ]R_|_ x M, and 
hence (|3.12p and the dominated convergence theorem implies that 

/■oo roo 

hm E(e-^i^'+^^2^') = / / 9oo(Ai, A2) P(y„,Xo)(dyo, dxo) = 9oo(Ai, A2) 

JO J-00 

for (Ai, A2) G X M. Then, using the continuity of the function ]R_|_ x M 9 (Ai, A2) 1— )■ 5oo(Ai, A2) 
(which will be checked later on), the continuity theorem yields (i). 

Next we turn to prove (j3.12p . Since > and 6 > 0, by ()3.9p . using also that vt{Xi, A2) ^ 
for all t ^ and (Ai, A2) G x M (see Step 1), we have 

lim {-yoVt{Xi, A2) + ixoe"^*A2) = 0. 



Recall that 

5t(Ai, A2) = -a / Ws(Ai, A2)ds + zmA2 







Since 6 > 0, we have limt_^oo ''—I — = ^, and since vt{Xi, A2) ^ for all t ^ 0, (Ai, A2) G M_|_ x M, 
(see Step 1), by dominated convergence theorem and (|3.9p . we get 



/•t />oo 

lim / ^^(Ai, A2) ds = / z;s(Ai,A2)ds. 
*^°°Jo Jo 



13 



Indeed, for all t ^ and s ^ 0, 

|?;^(Ai, A2)l[o,t](s)| Vs{Xi,X2), 

and, by ([^ . 

/ i;,(Ai,A2)ds ^ M(Ai,A2) / (1 + s) max(e-2^^ e"^^) ds 

JO JO 

= M(Ai,A2)maxQ + l,^ + ^) <oo. 

The continuity of the function M4. x M 9 (Ai, A2) 1— 500 (Ai, A2) can be checked as follows. It will 
follow if we prove that for all s ^ 0, the function Vg is continuous. Namely, if A^"^ = (A^"'\A2"^), 
n G N, is a sequence in x M, such that lim„_j.oo A^'") = A, where A G M+ x M, then 
lim„_j.oo ^^s(A("^) = Us(A) for all s ^ 0, and, by 



t;,(AW) = z;,(Ai"\A5")) ^ M(Ai"\A5"))(l + s)max(e-2^^e-'"'), n G N, s ^ 0. 
Since the sequence A^'"), n G N, is bounded (being convergent), we have 

supM(aS"\a^"^) < 00, 

nGN 

and using also that /o°°(l + s) max(e~^^*, e"'"*) ds < 00, the dominated convergence theorem implies 
that 

poo poo 

lim / vs{x[''\xi''^)ds= / ?;,(Ai,A2)ds, 
"^°°Jo Jo 

which shows the continuity of ^oo- Finally, we turn to prove that for all s ^ 0, the function 

R_l_ X M 3 (Ai,A2) I—)- fs(Ai,A2) is continuous. Note that the function Vg does not depend on 

the parameters a and m, since it is the unique solution of the differential equation (13. 3p . Let 

(Yt, Xt)t^o be an affine process satisfying the SDE (jl.ip with initial value {Yo,Xq) = {Yo,Xq) and 

with parameters a = m = and the given 6 > and ^ > 0. Then, by 



(3.13) e-^'^'+'^^^^Ps{{yo,xo),dO = exp {-yo^^.(Ai, A2) + ixoe-^^A2} 

for s G M+, (Ai,A2), {yo,xo) G M+ x M, where {Pt)t^o denotes the transition semigroup of the 
affine process {Yt, Xt)t^o. For all s G ]R_|_, the left-hand side of ()3.13p is continuous as a function of 
(Ai, A2) G M+ X M, since for ah (Ai, A2) G M+ x M and (/ii, /12) G M^, 

/ / e"(^i+'^i)«i+*(^2+'^^)«^P,((yo,xo),dO- / / e"^i5iW2p^((y^^^^)^dO 

Jo J -00 Jo J -00 

<; / / |e-'^^«i+*'^^«2-l|P,((yo,Xo),dO, 
JO J-00 

which tends to as (/ii,/i2) (0,0), by dominated convergence theorem. This implies that the 
right hand side of ()3.13p is also a continuous function (Ai,A2) G M+ x M. This readily yields the 
continuity of the function Vs for all fixed s G M+. 
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(ii): First we check that the one-dimensional distributions of (Yt, Xt)t^o are translation in- 
variant and have common distribution as (l"ooi^oo) has. Using (j3.2p . ()3.6p . the tower rule and the 
independence of {Yq,Xo) and {L,B), it is enough to check that for all t^O and (Ai, A2) G IR+ x M, 



E 



^exp|-'(;i(Ai,A2)yoo + ^*A2Xoo A2)|^ = exp |-a j Vs{Xi, X2) ds + i^X2^ ■ 



By ([321), (1121) and using also that vt{Xi, X2) ^ for all t^O and (Ai,A2) G M+ x M (see Step 1 
of the proof of part (i)), we have 

E (exp {-vt{Xi,X2)Y^ + ie-^*A2Xoo + gt{Xi, A2)}) 

= exp|-a j ■Us(i;i(Ai,A2),e"''*A2)ds-Fiye"^*A2 + 5t(Ai,A2)| 

= exp|-a^y" Vs{vt{Xi,X2),e"^^X2)ds + j ^^(Ai, A2) ds^ + i^A2| . 



ft 

-et ' 



Hence it remains to check that 

) /•oo 

Vs{Xi,X2)ds = / fs(^^i(Ai,A2),e^^*A2)ds-F / Vs{Xi,X2)ds 
Jo Jo 

for t ^ 0, (Ai, A2) G M+ X M, i.e., 

) /"OO 

VsiXuX2)ds= Vs{vtiXuX2),e-^'X2)ds, t ^ 0, (Ai, A2) G I 
Jo 

For this it is enough to check that 

?^s(^^t(Ai,A2),e~^*A2) = t;.+t(Ai,A2), s,t^O, (Ai,A2)gM+x 

or equivalently 

(3.14) vt{vs{XuX2),e-'>'X2)=Vs+t{XuX2), s, t ^ 0, (Ai, A2) G M+ x 

By (|3.3p . we have 



dv 



s+t 



(Ai,A2) = -bvs+tiXi,X2) - -{vs+t{Xi,X2)r + V2«(^+*)a2, t ^ 0, 



at a 2 

with initial condition Ws+o(Ai, A2) = ^^(Ai, A2). Note also that, again by (j3.3p . 

^(i;,(Ai,A2),e-^^A2) = -fe7;t(^;.(Ai, A2), e-^^A2) - -(^;t(^;s(Ai, A2), e-^^A2)r 
at a 

^l^_2ei^^_e.^^)2^ t;^0, 

with initial condition fo(ws(Ai, A2), e~^^A2) = Ws(Ai,A2). Hence, for all s ^ 0, the left and right 
sides of (j3.14p . as functions of t ^ 0, satisfy the differential equation ()3.3p with the initial value 
?;s(Ai,A2). Since (j3.3p has a unique solution for all non-negative initial values, we obtain ()3.14p . 

Finally, the strict stationarity (translation invariance of the finite dimensional distributions) of 
{Yt, Xt)t^o follows by the chain's rule for conditional expectations using also that it is a time homo- 
geneous Markov process. □ 
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4 Ergodicity 



Such as the existence of a unique stationary distribution, the question of ergodicity for an affine process 
is also in the focus of current investigations. 

Recently, Sandric |32j has proved ergodicity of so called stable-like processes using the same tech- 
nique that we applied. Further, the ergodicity of the so-called a-root process with a G (1,2] (see, 
the first SDE of (jl.ip ) and some statistical applications were given in Li and Ma |25| . 

The following result states the ergodicity of the affine diffusion process given by the SDE (|1.1|) 
with a = 2. 

4.1 Theorem. Let us consider the 2- dimensional affine diffusion model (jl.ip with a = 2, a > 0, 
6 > 0, m G M, 9 > 0, and with a (random) T^-measurahle initial value {Yq,Xq) such that 
P(yo ^ 0) = 1. Then for all Borel measurable functions / : — t- M such that E \ f{Yoo, Xoo)\ < oo, 
we have 

(4.1) ¥(^lim^^ £ fiYs,Xs)ds = EfiY^,X^)^ =1, 

where the distribution of (yoo^oo) is given by (j3.2p and (j3.3p with a = 2. 

Proof. We will use the notations of Meyn and Tweedie [27], [28]. Using Theorem 6.1 (so called 
Foster-Lyapunov criteria) in Meyn and Tweedie [28], it is enough to check that 

(a) {Yt,Xt)t^o is a right process (defined on page 38 in Sharpe [35]): 

(b) all compact sets are petite for some skeleton chain (skeleton chains and petite sets are defined on 
pages 491, 500 in Meyn and Tweedie [27], and page 550 in Meyn and Tweedie [26], respectively); 

(c) there exist c, d E M with c > such that the inequality 

{AnV){y,x) ^-cV{y,x) + d, {y,x)£On 
holds for all n G N, where 0„ := {{y,x) G M+ x M : ||(y,2;)|| < n} for each n G N, 
(4.2) v{y,x):={y-cif + {x-C2f, {y,x)£R+xR, 

with some appropriate ci,C2 G M, and An denotes the extended generator of the stopped 
process (Y'/"'^ Xj"^)f3.o given by 



(n) y{n). __ 



(YuXt), for t<Tn, 
(0,n), for t ^ r„. 



where the stopping time r„ is defined by T„ := inf{t G M+ : (Yt,Xt) G (M+ x M) \ 0„}. 

Indeed, then Theorem 6.1 in Meyn and Tweedie [28] yields the exponential ergodicity of the process 
{Yt,Xt)t£M.+ , namely, there exist /3 > and G M+ such that 



sup 



E {g{Yt, Xt) I (Fo, ^o) = (yo, xq)) - E(g(yoo, ^oo)) ^ B{V{yo,xo) + l)e 
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for all t G M+ and {i/q, xq) £ x M, where the supremum is running for Borel measurable functions 
g : M_|. X M — >■ M. According to the discussion after Proposition 2.5 in Bhattacharya [5], this implies (i). 
Here we also point out that, due to Bhattacharya [5], we do not have to assume that ^{Yq > 0) = 1 
in order to prove (i), we will need this additional assumption only for part (iii). 

To prove (a), it is enough to show that the process (Yt, Xt)t£R_^ is a (weak) Feller (see Meyn and 
Tweedy [271 Section 3.1]), strong Markov process with continuous sample paths, see, e.g., Meyn and 
Tweedy [271 page 498]. According to Proposition 8.2 (or Theorem 2.7) in DufHe et al. |10| . the process 
{Yt,Xt)t^o is a Feller Markov process. Since (Yt,Xt)t^o has continuous sample paths almost surely 
(especially, it is cadlag), it is automatically a strong Markov process, see, e.g., Theorem 1 on page 56 
in Chung [8]. 

To prove (b), in view of Proposition 6.2.8 in Meyn and Tweedy [29], it is sufficient to show that 
the skeleton chain (Yn, Xn)nez+ is irreducible with respect to the Lebesgue measure on x R, 
and admits the Feller property. The skeleton chain {Yn, Xn)nez+ admits the Feller property, since 
the process (Yt, Xt)t^o is a Feller process. In order to check irreducibility of the skeleton chain 
ez+ with respect to the Lebesgue measure on M+ x M (see, e.g., Meyn and Tweedy [281 page 
520]), it is enough to prove that the conditional distribution of {Yi,Xi) given {Yq,Xq) is absolutely 
continuous (with respect to the Lebesgue measure on M+ x R) with a conditional density function 
/(yi,Xi)|(yo,Xo) : X R2 ^ R+ such that | (yo,Xo)(yi ^ I yo, a^o) > for all {y,x,yo,xo) G 

(0, oo) X R X R+ X R. Indeed, the Lebesgue measure on R_|- x M is cr-finite, and if i? is a Borel 
set in R+ x M with positive Lebesgue measure, then 

E(f2lBiYn,Xn) (Yo,Xo) = {yo,xo)] ^P((yi,Xi)GS|(yo,^o) = (yo,xo)) 

= jj f{Y^,Xi)\{Yo,Xo)iy^^\yo,xo)dydx > 

B 

for all (yo; a^o) £ 1^+ X I^- The existence of f(Yi,Xi)\{Yo,Xo) with the required property can be checked 
as follows. By Theorem 12.21 we have 

Yi = e-'' (^yo + aj^ e^" du + e^'^y^dL^^ , 

Xi = e-^ (^xo + m e'^^'du + J^ e^"v^dS„ 

provided that {Yq,Xq) = {yo,xo), {yo,xo) G R+ x R. Using that a two-dimensional random vector 
is absolutely continuous if and only if + v is absolutely continuous for all invertable matrix 
V G M^^^ and for all vector t; G R^, it is enough to check that the random vector 

(4.3) (^j\'>^V^dLu, j^'^'y^dBu 

is absolutely continuous with respect to the Lebesgue measure on R^ having strictly positive density 
function on R^. Using the independence of the standard Wiener processes L and B, we have for 
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all X, y G 



E P 



E 



0^t)t&[o,i] 



y i-x 



oo J — oo 

y fx 



piz; I e^^^'Yudu ] p { w; I e'^^'Y^du ) dzdw 



oo J — oo 



Elplz;/ e2^"y„ du ) p ( / e^^^Y^ d-u ) ) dz 



where 



1 



> 0. 



Consequently, the random vector (j4.3p is absolutely continuous having strictly positive density 



9 {w, v)^W.(^ l^- e^'^Yu dn^ p e^^^y. 



du 



Here we call the attention that due to the assumption a > 0, 



e^'^Yu du>0 



e^^^y^^dn > ) = 1. 



We remark that Ben Alaya and Kebaier [H Proposition 2] derived an explicit expression for the density 
function of Yt, t > 0. 

To prove (c), first we note that, since the sample paths of (Y,X) are almost surely continuous, 
for each n G N, the extended generator has the form 

{■^nf){y,x) = ^yfli{y,x) + ^y/a'sl^ja;) + (a - by)f[{y,x) + (m - 9x)f[{y,x) 

for all {y,x) G On and / G C'^{R+ x M,]R), see, e.g., page 538 in Meyn and Tweedie [28]. We also 
note that, by Duffie et al. [IHl Theorem 2.7], for functions / G C^{R+ x ^^^)^ AJ = Af on 0„, 
where A denotes the (non-extended) generator of the process {Yt,Xt)t&R^- For the function V 



defined in we have V G C\^+ x M,R) and 

= 2(y-ci), Vi{y , x) = 2{x - ci) , Vi[^{y , x) = V^^iy , x) = 2 

for (y, x) G M+ x M, and hence for all n G N and < c < oo, 

{AnV){y,x) + cV{y,x) 
= 2y+{a- by)2{y - ci) + (m - 6'x)2(x - C2) + c{y - cif + c{x - 02)^ 
= (c - 2b)y'^ + 2(1 + a + bci - cci)y + ccj - 2aci + 2(m - 9x){x - C2) + c(x - 02)^ 
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for all (y, x) G 0„. Let us choose c > in such a way that 2b and let C2 := Then 
{AnV){y,x)+cV{y,x) 

= (c - 26)y^ + 2(1 + a + 6ci - cci)?/ + ccf - 2aci -26'(x-— ) +c(x- — 



1 + a + 6ci — cci \ (1 + a + 6ci — cci) ^ 2 



= (c-25)^, + ^3^6 J c-26 

Hence, 

{AnV\y,x)^-cV{y,x)^d, {y,x)£On, n G N, 
with c G (0,2min(6, 6*)) and 

(l + a + 6ci -cci)2 2 o 

d := — h cci - 2aci. 

c — 2b 

Note that the above argument also shows that in the definition ()4.2p of the function V the constant 
ci can be an arbitrary real number. □ 

In the next theorem we collected several facts about the random variable ( Yqo > ^00 ) given by 
(13:2]) and (1331) with a = 2. 



4.2 Theorem. Let us consider the 2- dimensional affine diffusion model (jl.ip mt/i a = 2, a > 0, 
6 > 0, m G M, > 0, and with a (random) To-measurable initial value (Yq^Xq) such that 
P(yo ^ 0) = 1. Then the random variable (1^00) -'^oo) given by (|3.2p and ()3.3p mt/i a = 2 is 
absolutely continuous, the Laplace transform of Y^o takes the form 



\ \ —2a 

(4.4) E(e-^i^-)= (l + ^J , A1GM+, 

yielding that has Gamma distribution with parameters 2a and 2b. Further, all the (mixed) 

moments of (YooiXoo) of any order are finite, i.e., we have ^{Y^\XooY') < °o for all n,p G 
and especially, 

a iTi 

E(yoo) = nxoo) = J, 

E(yoo^i) = (6 + 2g)262g2 (^(«^ + + ^) + 2"^'K2e + &))• 

Proof. First we show that the random variable ((Yoo-,Xoo) is absolutely continuous. Let us suppose 
that the random initial value {Yq,Xq) has the same distribution as that of {Yoo,Xoo)- Then, by 
part (ii) of Theorem 13.11 the process (Yt,Xt)t^Q is strictly stationary. Hence it is enough to prove 
that {Yi,Xi) is absolutely continuous. This can be done similarly to the proof of condition (b) in 
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the proof of Theorem 14.11 we only sketch a proof. Namely, by Theorem 12.21 we have 



¥1=6-" [Yo + a I e^" du + 



Xi = e-" [Xo + 7n I e^" du 



Using that if ^ and rj are independent real- valued random variables such that one of them is absolutely 
continuous, then their sum S, + rj is absolutely continuous (regardless of the other distribution), it is 
enough to check that the random vector 



e^''^/Y,dB, 



is absolutely continuous with respect to the Lebesgue measure on M?. Its proof goes along the very 
same lines as in the proof of condition (b) in the proof of Theorem 14.11 

Now we turn to prove that Yoq has Gamma distribution with parameters 2a and 2b. By (|3.2p . 



Efe" 



exp < —a 



i;s(Ai,0) ds 



Ai G 



where vt{Xi,0), t^O, is the unique solution of the differential equation 
(4.5) 



|^(Ai,0) = -bvt{Xi,0) - UvtiXuO))^ t ^ 0, 



[vo(Ai,0) = Ai. 

Note that the formula above for E(e~'^i^°°) is a special case of formula (3.31) in Li [2l]. The 
differential equation (j4.5p is of Bernoulli type, and hence for Ai > 0, with the transformation 
iii(Ai,0) := Ut(Ai,0)~^, t ^ 0, we have the inhomogeneous linear differential equation with constant 
coefficients: 



(4.6) 

One can check that 



■^(Ai,0) = 6nt(Ai,0) + i, t^O, 



[7xo(Ai,0) = Ar^ 
..(A„0)=..(A„0)- = ((i- + l)e--l 



-1 



Ai > 0, t ^ 0. 



Hence, 



exp < 



exp < 



^00 J 

io 7X711 



Ai + 26 



+ 



phs _ J_ 
^ 2b 



■ds 



1 + 



26 



exp < 



dz 



dz 



-2a lim ( — In z + In 

2—^00 
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which yields (j4.4p . 

Next we give another argument for deriving (j4.4p . Since in case of a = 2, the process (Yi)t>o is a 
continuous state branching process with branching mechanism hz + z^ /2, z ^ 0, and with immigration 
mechanism az, z ^ 0, by the proof of Theorem 3.20 in Li |24j . we have 



/ 

Jo 



v.iX„0)as = l ^J^d. = 2aln^l + -), A, E 



'0 ^ 
which yields ()4.4p . 

Now we prove that all the mixed moments of {Yao, X^o) are finite. First we show that 



(4.7) 



/ E{Y,^'XlP)du < oo for aU t^O and n,p G 



whenever all the mixed moments of the initial value (Yq,Xq) are finite and P(Yo > 0) = 1. Due to 
the inequality between two power means, (a + 6 + c)^^ ^ K{o?p + iP'^ + c^^), a,b,c £ M, with some 
constant K > 0, and hence, by (12. 2p . to prove ()4.7p it is enough to check that 

(4.8) E (y^ ^e-^f^^Xo^P + (^m e'^^""^) dt;^ + (^J^" e-^^""^') v^dS,,^ "^Au < oo 

for all t^O and n,pGZ+. Since for all ti€ [0, the conditional distribution of JJ* e~^^"~^^-v/T^d-B„ 
with respect to the u-algebra generated by (ys)sg[o,t] is a normal distribution with mean and with 
variance Jq" e"^^*-""^)!^, dz;, for proving ()4.8p it is enough to check that for all t ^ and n,p € Z+, 

/* E(e-2P^"y„"Xo2^) du < oo, /* E(y„") dn < oo, 
7o Jo 



and 

ft 



/ e-2^("-^)K,df 1 |d^x<oo. 



Using that > 0, "^{Yt ^ 0, t £ IR+) = 1, by Cauchy-Schwarz inequality, for all u £ [0,t], we have 



1/2 



ajid 



1/2 



1/2 



1/2 



mul) " ' • • • ^" E (n, • • • i;, J dz;i • • • dr;2p) 

^ ( E (e)) (e (Y^'r) • • • E )) '^^''^ . . . d.2p) , 

where the last inequality follows by the multivariate version of Holder's inequality. Since E{X^P) < oo, 
this shows that in order to prove (|4.7p it is enough to check that for all k £ Z+ and t £ M+, 

sup E{Y^) < oo. 

ue[o,t] 
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According to Proposition 3 in Ben Alaya and Kebaier [3], 

(4.9) supE(y/ I To = yo) < oo for all yo > and k G Z+. 

Next we show that for all s ^ 0, yo > 0, and k G 

= yo) is a polynomial of yo of degree k. 
By (HSl), for all s 5? 0, Ai ^ 0, and (yo,2;o) G M+ x M, 

/•oo POO r 1 

(4.10) / / e-^i«ip,((yo,xo),de)=exp<^-yo?;,(Ai,0)-a / i;„(Ai,0)dnL 

Since, by (fOj) . E(y^^ | Iq = 2/o) < oo for all yo > and k e and, by (j33]) . the function 

i;s(Ai, 0), s ^ 0, is continuously differentiable of infinitely many times, one can differentiate both sides 
of ()4.10p with respect to Ai for k times. Since ^^(0, 0) = for all s ^ (which is a consequence 
of the uniqueness of the solution of the differential equation ()3.3p ). we have 

E{Y^^ I yo = yo) = (-1)'=^ exp |-yo^;.(Ai, 0) - a Vu{Xi,0) du 

is a polynomial of yo of degree fc, and the coefficients of this polynomial are continuous functions 
of s. Since a continuous function on a compact set is bounded, we get that sup„g[o,t] '^^u I ^ = yo) 
can be bounded above by a polynomial of yo having degree /c, say Qkivo) (this polynomial depends 
also on t). Hence, by the law of total expectation, 

sup E(y„'=) = sup / E(y„'=|yo = yo)Py„(dyo) ^ / Qfe(yo) Pi^, (dyo) = E(Qfc(yo)) < oo, 

ue[o,t] u&[o,t] Jo Jo 

where the last inequality follows by the assumption that all the (mixed) moments of (yo,Xo) of any 
order are finite. 

We note that for proving the finiteness of K{Y^Xu^) du we could have used part (i) of Theorem 
2.16 in Duffie et al. [10]. This way of proving is somewhat complicated that's why we decided to find 
another way presented above. 

For all n,p £ using the independence of L and B, by Ito's formula, we have 

d{Yt"Xf) = nYt^~'Xf{{a - bVt) dt + VVtdLt) +pY;'Xf-\{m - OXt) dt+^tdBt) 

+ ^^^^Yr'xfY, dt + ^-^^Yrxr'Y, dt 

nYr\a - hY,)Xl+pYr{m - eX,)Xr' + "^^^^Yr'Xf 



2 

for t ^ 0. Writing the SDE above in an integrated form and taking expectations of both of its sides, 
we have 

mrxf) - E(yo"Xo^) = f lanEiYr'XZ) - &nE(yj^XS) +pmE{Y:XZ~') -p^E(yj'X^) 







du, t > 0, 
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where we used that 



/* Vr^l-'Xl dL^) and { f Y^+^'XP-' dB,, 

Jo / t^o \Jo 



t>0 



are continuous square integrable martingales due to (|4.7p . see, e.g., Ikeda and Watanabe [T71 page 55]. 
Introduce the functions fn,pit) :=E{Y^^Xf), t G IR+, for n,pGZ+. Then we have 

/^pW = -ibn+pe)fnAt) + (an + ifclll^ + pmU,p-iit) + ^^^^/„+i,p_2(t) 

for t e M+, where fk,e{t) := if k,£ e Z with < or £ < 0. Hence for all M G N, the 
functions fn,p, n,p £ with n + p ^ M satisfy a homogeneous linear system of differential 
equations with constant coefficients. For example, if M = 2 then 





























/{,o(*) 
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-b 
















/i,o(t) 


Alw 




m 





-9 













hi{t) 


fioit) 







2a + 1 





-26 










/2,0(t) 


fllit) 







771 


a 





-b-e 







flAt) 


/o,2(i). 
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2m 








-29 




/0,2(t) 



t £ 



Thus, for all n,pGZ+, the function is a linear combination of the functions Q-{^^+^S)t ^ t G ]R_|_, 
k,i € with A; + £ ^ n + p, since the eigenvalues of the coefficient matrix of the above mentioned 
system of differential equations are —{kb + i9), k,i G Z_|_ with k + I ^ M. Consequently, for 
all n,p G Z+, the function fn^p is bounded on R_(_ and the limit limt^ao fn.p{t) exists and 
finite. By the moment convergence theorem (see, e.g., Stroock [3U Lemma 2.2.1]), lim^^i^oo fn,p{t) = 
limt_>oo ) = 'KiY^X'^)^ n,p G Z+. Indeed, by Theorem 13.11 and the continuous mapping 

theorem, Yfl^Xf — > Y^X^ as t — t- oo, and the uniform integrability of {Y^"Xf : t G follows 
from the boundedness of the function f2n.2p on M+. Hence we conclude that all the mixed moments 
of {Y^,X oo) are finite. 

Finally we calculate the moments listed in the theorem. Suppose that the initial value (Yq,Xq) 
has the same distribution as {Yoo,Xoo)- Then, by Theorem 13.11 the process {Yt,Xt)t^o is strictly 
stationary, and hence, fn,p{t) = lE(y^X^), for all t G M4. and n^p G Z_|_. The above system of 
differential equations for the functions /„,^p, n,p G Z4,, yields 



1 



bn + p9 



(4.11) 



an + ] E(y^-iX^) + mpE{Y^XP^'^ 



for all n,p G Z-|_. By (j4.1ip . one can calculate the moments listed in the theorem. 

Finally, we note that for calculating the moments E(Y^X^) < 00, n,p £ Z_|_, we could have used 
formula (4.4) in Filipovic et al. [12] which gives a formal representation of the polynomial moments 
of (Yt,Xt), t G M+. The idea behind this formal representation is that the infinitesimal generator 
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of the affine process (Y, X) formally maps the finite-dimensional linear space of all polynomials in 
(y, x) £ M+ X M of degree less than or equal to k into itself, where G N. We also remark that the 
moments of Yoo could have been calculated directly using that Y^q has Gamma distribution with 
parameters 2a and 2b. □ 
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